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Abstract 



A search for groups of collimated muons is performed using a data sample collected 
by the CMS experiment at the LHC, at a centre-of-mass energy of 7 TeV, and corre- 
sponding to an integrated luminosity of 35 pb^^. The analysis searches for production 
of new low-mass states decaying into pairs of muons and is designed to achieve high 
sensitivity to a broad range of models predicting leptonic jet signatures. With no ex- 
cess observed over the background expectation, upper limits on the production cross 
section times branching fraction times acceptance are set, ranging from 0.1 to 0.5 pb 
at the 95% CL depending on event topology. In addition, the results are interpreted 
in several benchmark models in the context of supersymmetry with a new light dark 
sector exploring previously inaccessible parameter space. 
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1 Introduction 

Recent astrophysical observations of an excess of high-energy positrons in the cosmic-ray spec- 
trum 111 have motivated the rise of new physics scenarios IH suggesting that this excess may 
be associated with annihilations of dark matter particles in the galactic halo. These models may 
also accommodate the observed discrepancies in direct searches for dark matter [3 -5J. One re- 
alization of such models assumes an extra U{1) gauge symmetry with weak coupling to the 
standard model (SM). The U{1) symmetry is broken, leading to a light vector boson (with mass 
m ~ 0(1 GeV/c^)), referred to as a "hidden" or "dark" photon 7dark/ which can have a small 
kinetic mixing with the SM photon providing a portal for the hidden-sector photon to decay 
into leptons and, if kinematically allowed, hadrons. More complex models can lead to a whole 
hierarchy of hidden-sector states or can have dark photons preferentially couple to leptons ||6|. 
Hidden sectors can be realized naturally in supersymmetric (SUSY) models where coupling of 
the dark sector to the SUSY sector can be enhanced. At the Large Hadron Collider (LHC), if 
SUSY exists and is kinematically accessible, these models predict production of dark photons 
as part of the SUSY cascades. The new light hidden states may be produced in decays of the 
lightest supersymmetric particle (sparticle) of the visible sector. In the extensions of the Min- 
imal Supersymmetric Standard Model (MSSM) with the lightest neutralino being the lightest 
SUSY particle (LSP) in the visible sector, the MSSM LSP can decay to the light hidden-sector 
particles and sparticles, while a heavy dark fermion provides a cold dark matter candidate. Al- 
ternatively, the lightest MSSM neutralino can decay to light hidden particles and a heavier dark 
neutralino [7J, which becomes a cold dark matter candidate. Because the MSSM LSP in these 
models is unstable and is not subject to the constraints on a cold dark matter candidate, there 
are scenarios [8] where the LSP is a squark decaying into a quark and the light hidden-sector 
states. Depending on the complexity of the light dark sector, at the LHC one may expect either 
a single dark photon at the end of each SUSY cascade or a whole cascade of hidden-sector-state 
decays with emission of multiple dark photons. Subsequent decays of new states into leptons 
lead to appearance of energetic collimated groups of leptons, the characteristic "leptonic jet" 
signature |9|. While the spectrum of such leptonic jets depends on masses of sparticles that are 
part of the SUSY cascade in which they are produced, typical scenarios accessible at the LHC 
predict leptonic jets with transverse momenta ranging from tens to several hundreds of GeV/ c. 

Previous searches for low-mass dilepton resonances have been performed at the Tevatron llTOl 
[TT| , as well as Belle |(T2|, BaBar |13|, and LEP tMH , and revealed no signals of new physics. 
Due to the large cross section of SUSY production via strong interaction, the LHC may provide 
access to these new states with early data motivating the search for anomalous production of 
collimated groups of leptons. 

This paper describes a topology-based search for new light resonances decaying to pairs of 
muons using data collected by the Compact Muon Solenoid (CMS) experiment during the 2010 
LHC data-taking period and corresponding to an integrated luminosity of 35.0 ± 1.4 pb^^. This 
analysis searches for events with one or more muon pairs that are consistent with being pro- 
duced in the decays of the same particle type. Assuming on-shell production of at least a frac- 
tion of these new bosons per event, the new physics would manifest itself as an enhancement 
in the production rate of muon pairs consistent with a certain common mass. The results of the 
search are presented in a model-independent fashion, limits are also set on specific benchmark 
scenarios LZ^ in the context of SUSY models. 
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2 Detector, Dataset, and Trigger 

The CMS detector is a general-purpose apparatus providing excellent momentum and direc- 
tion measurement of particles produced in pp collisions at the LHC. The central feature of 
CMS is a superconducting solenoid of 6 m internal diameter. Within the volume of a 3.8 T 
magnetic field are the silicon pixel and strip tracker, the crystal electromagnetic calorimeter 
and the brass/ scintillator hadron calorimeter. Muons are measured in gas-ionization detectors 
embedded in the steel return yoke. In addition to the barrel and endcap detectors, CMS has 
extensive forward calorimetry. CMS uses a right-handed coordinate system, with the origin 
at the nominal interaction point, the x axis pointing to the centre of the LHC ring, the y axis 
pointing up (perpendicular to the LHC plane), and the z axis along the counterclockwise beam 
direction. The polar angle is measured from the positive z axis and the azimuthal angle (p is 
measured in the x-y plane. The pseudorapidity tj = — ln[tan(0/2)] is frequently used instead 
of the polar angle 0. Here we only briefly describe the components of the CMS directly relevant 
to this analysis; the full details of the detector, its subsystems, and performance are described 
elsewhere |(T5|. 

The inner tracker measures charged particles within the pseudorapidity range \tj\ < 2.5. It 
consists of 1440 silicon pixel and 15 148 silicon strip detector modules and is located inside 
the superconducting solenoid. It provides an impact parameter resolution of ~ 15 /^m and a 
transverse momentum (pj) resolution of about 1.5% for 100 GeV/c particles. The muons are 
measured in the pseudorapidity window |?/| < 2.4, with detection planes made of three tech- 
nologies: Drift Tubes, Cathode Strip Chambers, and Resistive Plate Chambers. Matching the 
muons to the tracks measured in the silicon tracker results in a transverse momentum resolu- 
tion between 1 and 5 %, for pj values up to 1 TeV/ c. 

Because of the high rate of the collisions, the CMS uses a two-level dedicated trigger system. 
The first level (LI) of the CMS trigger, composed of custom hardware processors, is designed 
to select, in less than 1 }is, the most interesting events, using information from the calorimeters 
and muon detectors. The High Level Trigger (HLT) processor farm, running a simplified and 
highly optimized version of the CMS offline reconstruction, is designed to decrease further the 
event rate down to a maximum of 300 Hz, before data storage. 

The data used in this analysis have been collected with the inclusive muon triggers with the 
lowest available pj threshold. At level 1, the data are selected using muon candidates recon- 
structed by the LI muon hardware, followed by the confirmation at the HLT, where muons 
are reconstructed by matching tracks reconstructed in the muon system with the tracks recon- 
structed in the silicon tracker detectors to refine the muon pj measurement. Because the trigger 
configuration was changing during the data taking period, there are three distinct parts of the 
dataset where the triggers used had transverse momentum thresholds of 9, 11, and 15 GeV/c 
at the HLT. In all cases, level 1 muon thresholds were low enough to ensure that the HLT 
thresholds are at the plateau of the LI muon efficiency. To make the selected data uniform, we 
additionally require offline events to contain at least one trigger candidate with pj > 15 GeV/ c 
as measured online, such that the final dataset is the same as it would have been if collected 
using a single inclusive muon trigger with pj > 15 GeV/ c. 

3 Offline Selections and Analysis of the Data 

In this analysis we search for evidence of new light bosons decaying into pairs of muons. The 
new particles can be isolated or produced in groups, coming from cascade decays in the hid- 
den sector ending with several instances of the lowest-state particle. In addition to dimuon 
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decays, the new particle decay channels can include pairs of electrons and perhaps hadrons. 
In many scenarios, multiple instances of such a new boson can be produced per event. This 
analysis therefore searches for events with one or more muon pairs that are consistent with 
being produced in the decays of the same particle type. Assuming on-shell production of at 
least a fraction of these new bosons per event, the new physics would manifest itself as an 
enhancement in the production rate of muon pairs consistent with a certain common mass. 

This analysis aims to reach two goals. First, to achieve high sensitivity for a representative 
range of specific new physics scenarios leading to characteristic muon-jet signatures and, sec- 
ond, to present results that would allow future interpretation in the context of other models of 
new physics yielding leptonic jet signatures. Essentially all classes of models of new physics 
relevant to this analysis lead to the production of events of several different topologies in terms 
of the number of coUimated muon jets and multiplicities of muons within each jet. Because 
these topologies have different sources and levels of SM backgrounds, we categorize events 
using topologies with different expected signal-to-backgroimd ratios to maximize the overall 
sensitivity of the analysis. 

3.1 Data Selections and Reconstruction 

Offline, events are required to have at least one primary vertex reconstructed in the luminous 
region along the beamline to minimize background events not originating from collisions. 
Selected events are further required to have at least one high-quality muon candidate, with 
> 15 GeV/c matching the muon selected by the online trigger and within |//| < 0.9 re- 
constructed using an "inside-out" algorithm, thus ensuring high efficiency in the environment 
with multiple nearby muons. This algorithm extrapolates silicon tracker tracks into the muon 
system and attaches to them individual tracklets (stubs) reconstructed in muon chambers. Any 
stub in the muon system can only be associated with one extrapolated tracker track most com- 
patible with the stub. The reconstructed muon candidate is required to have stubs in at least 
two out of four muon stations it crosses. To be classified as high-quality, muon candidates are 
required to have at least eight hits in the silicon tracker. The requirement of |//| < 0.9 ensures 
high and well understood trigger efficiency insensitive to the presence of muon hits from other 
nearby muons expected in the signatures with collimated muons. It avoids the endcap region 
where the trigger efficiency can be substantially diminished in the presence of multiple closely 
spaced muons because of the features of the trigger electronics setup. Additional muon candi- 
dates are required to have pj > 5 GeV/ c, to be contained within | j/| < 2.4, and to satisfy the 
same quality requirements. 

Muon jets are reconstructed by iteratively clustering muon candidates starting with the highest 
Pt- Each additional muon is added to the jet if the invariant mass of this muon and any op- 
positely charged muon already assigned to the jet satisfies < 9 GeV/ and is compatible 
with originating from the same vertex (confidence level of the vertex fit > 1%). The cluster- 
ing procedure always converges and is independent of the order in which muons are added 
to muon jets. The choice of mf^f^ < 9 GeV/ ensures that muons originating from the same b 
quark are always clustered into the same muon jet and most muons originating from different 
b quarks are clustered into different muon jets. It is also appropriate for topologies predicted 
by most relevant models of new physics, as the typical masses of the heavier hidden states 
originating the cascades are of the order of a few GeV/ c^. Note that the efficiency of this clus- 
tering algorithm does not depend on the boost of the muon jet, thereby reducing the sensitivity 
of the analysis to those details of the kinematics that can differ from one model to another. As 
a consequence of the clustering algorithm, each muon jet must contain at least one muon of 
each charge, but can contain arbitrarily many muons. Within each muon jet, we identify "fun- 
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damental dimuons," the pairs of oppositely charged muons that are most likely to have arisen 
from individual dark photon decays. Since all dark photons in the event have the same mass, 
the assignment is performed by selecting a combination of pairs that yields minimal difference 
in dimuon mass among the pairs. 

3.2 Event Categorization 

All events are categorized according to the number of muon jets N and the number of muon 
candidates in the muon jet, thus forming topologies denoted as i^^^..,,^. No isolation 
requirements are imposed on any of the muon candidates. Muons that do not belong to any 
muon jet (which may arise from SUSY cascades, rather than dark photon decays) are neither 
used to identify nor to reject signal events. The event selection is also insensitive to the presence 
of hadronic jets or missing transverse energy in the event. 

In topologies with multiple dimuon candidates originating from the same particle type, the 
reconstructed masses of all fundamental dimuons would be consistent with each other within 
detector resolution for signal events, but not necessarily for backgrounds. Therefore, we build a 
K-dimensional distribution of reconstructed dimuon masses nii, niK for each topology, where 
K is the number of reconstructed dimuons per event. The signal of new physics would appear 
as an enhancement of events at a point near the K-dimensional diagonal with nii ~ ... ~ thk ~ 
mo, where mq is the mass of the new particle. While the distributions of backgroimd events 
are not smooth because of the low-mass resonances, they extend beyond the diagonal in a 
known way. The only exception is the R\ topology with exactly one fundamental dimuon per 
event: the signal would appear as a narrow peak in the ID distribution of dimuon mass m. 
For topologies with K dimuon candidates, we define the signal region as a "corridor" near the 
diagonal in the X-dimensional space of width 5 x (7(m„^,), where C7(m^^,) = 0.026 GeV/c^ + 
0.0065 m^^ for barrel {\t]^^^\ < 0.9, where t]>^^ is the psedorapidity of the dimuon momentum 
defined as a vector sum of momenta of the muons in the dimuon) and 0.026 GeV/c^ + 0.013 
for endcap (0.9 < \t]^^\ < 2.4). The parameterization for cr{mfip) was derived from studies of 
J/i/j, xp', (p, and p/ CO resonances, as well as high-pj Monte Carlo simulations, and corresponds to 
the resolution expected of hypothetical dimuons with « 300 GeV/ c in the barrel region and 
Pt ~ 150 GeV/ c in the endcap. Better mass resolution for lower momentum dimuons makes 
this definition of the signal corridor conservative for < 300 (150) GeV/ c in the barrel 
(endcap). Dimuon momentum spectum predicted in typical models rarely extends beyond 300- 
400 GeV/ c making this choice acceptable for the entire range of expected dimuon momenta. 
After the shape of the backgroimd events distribution in the K-dimensional space is measured, 
the data in the off-diagonal part can be used to obtain the background normalization, which 
can then be used to fit the data in the near-diagonal region for signal plus backgroimd. 

If an enhancement were to be observed in the diagonal regions, one could further construct 

the invariant mass of combinations of dimuons in the same muon jet to search for possible 
structure, e.g. a process of the type a2 — )■ flifli — )■ would lead to an enhancement 

in the invariant mass of pairs of dimuons consistent with m(fl2)- To maintain analysis sensi- 
tivity for models with m(fl2) < 2m(fli) where one of both fli bosons are produced off-shell, 
events with high multiplicity muon jets, e.g. topology R\, are to be analyzed for evidence of an 
enhancement in the "quadmuon" invariant mass distribution. 

3.3 Background Estimation 

The SM background contributions to event topologies jR^..„j^ vary significantly in terms of 
the overall rates and the composition of contributing processes. The single-dimuon topology 
R\ suffers from a particularly high rate of the SM backgrounds because of bb and Drell-Yan 
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Figure 1: Dimuon invariant mass distributions in the background-enriched samples for differ- 
ent topologies. The data are overlaid with parameterized functions, fitted to the data, and used 
to construct mass-shape templates for the distributions of background events in signal regions. 



processes. Without additional selections, the SM backgrounds would be too large to maintain 
sensitivity to signals with picobarn-scale cross sections. To reduce the SM backgrounds, events 
in the R\ topology are additionally required to have the muon jet transverse momentum satisfy 
Pj^' > 80 GeV/ c. This requirement dramatically improves the sensitivity of this topology to 
new physics signals that predict highly boosted muon jets. At the same time this requirement 
reduces acceptance for signal events containing only one dimuon per event, particularly for 
models, in which the muon jets have lower boosts. For such models the sensitivity is driven 
by events with topologies containing two or more dimuons per event, for which no high mo- 
mentum requirement is imposed. Data reconstructed in topology R\ with lower momentum 
dimuon candidates are used for background studies and validation of the background estima- 
tion techniques. The other special case is topology Kg, in which the signal content is expected to 
be low and the background rate is substantial, dominated by events with two genuine muons 
from b decays and a non-muon track misidentified as a muon. We do not use data in this topol- 
ogy to search for signal; these events are instead used as a control region to test background 
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parameterizations. 

The background rates in the selected topologies are expected to be low with the exception of 

topology R\, where backgrounds remain non-negligible even after the p^'' > 80GeV/c require- 
ment. However, because the search for new resonances is performed in small windows in the 
invariant mass distribution of muon pairs, the rate of remaining background in each window 
is comparable to the rate of the signal being sought. For topology R2, the main SM background 
contributions are due to bb production with one of the b quarks undergoing a double semilep- 
tonic decay, low-mass resonance production (prompt or from heavy-flavour decays), low-mass 
Drell-Yan production, and occasional muon misidentifications due to decays-in-flight, either 
alone or in combination with a muon from a heavy-flavour decay. For topology R22, the SM 
backgrounds are dominated by bb production with both b jets undergoing double semileptonic 
decays or fragmenting into low-mass resonances decaying to pairs of muons. Background 
events with muon jets consisting of multiple muon candidates (the quadmuon topology R\ 
and the higher-order topologies) typically originate from events with two muons from a b jet 
and the other muons from either decays-in-flight, punch-through, or muon misidentifications 
where some of the segments from true muons are matched to the non-muon tracks. The SM 
content of the higher order regions is due to rare combinations of the mechanisms discussed 
above and is extremely low. 

To account for background contributions, we construct "templates" (one for each topology) 
modeling the distribution of reconstructed muon pair masses in background events. With the 
exception of the single dimuon topology R\, the templates are multidimensional distributions 
in the (mi, mjc) space of reconstructed muon pair masses, where K is the number of dimuons 
characteristic of a given topology. For each category, we define one or more backgroimd- 
enriched regions used to construct the template. In addition, we define a control region for 
validating the template using events with properties closely resembling those of the final events 
or using the off-diagonal side-band of the final (mi, m^) distribution. 

While the templates were derived directly from data, we use simulation to determine the com- 
position of the backgrounds as well as momentum evolution of certain parameters, e.g. the 
dimuon mass resolution and shape of the invariant mass distributions. To ensure that simu- 
lation is reliable in the phase space characteristic of this analysis, a series of detailed studies 
have been performed. First, the single-dimuon dataset with < 80 GeV/ c (low momentum 
part of topology R]) was divided in subsets, each one dominated by only one of the contribut- 
ing background processes to measure rates, shapes, kinematic distributions, tracking related 
parameters and resolutions, and mass resolutions of low-mass resonances as a function of the 
boost. These measurements were compared to simulation predictions and showed very good 
agreement except for a few known and well-understood shortcomings of the simulated sam- 
ples available (lack of very low-mass Drell-Yan events, missing production modes and /or de- 
cay channels for some resonances). In the following, we describe the details of the procedure to 
construct the background templates for topologies R\, R22 and R\, which have the highest back- 
ground content. Templates for higher order topologies are constructed in a similar fashion, but 
ultimately have not been used explicitly in the final fit because no data events were observed in 
topologies with more than four muons (the Bayesian limit is insensitive to the rate of expected 
background events when no data events are observed). These topologies still contribute to the 
analysis sensitivity as benchmark models predict a non-negligible fraction of signal acceptance 
to fall into the higher order topologies. 

To model the shape of the invariant mass distribution for the single dimuon region R^, we 
define two sub-regions with 40 < < 60 GeV/ c (background-enriched region) and 60 < 
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Figure 2: (a): Data in the single-dimuon category {R\) control region 60 < p'^'' < 80 GeV/c 
overlaid with the background prediction obtained from the background-enriched region 40 < 
Pj^' < 60 GeV/ c, fitted for overall normalization only, (b): The invariant mass of all dimuons in 
the off-diagonal region for events in the two-dimuon category (R22' riote that there are two en- 
tries per event), compared with the prediction obtained from the full 2D background template, 
fitted for overall normalization only, (c): The invariant mass distribution of all "dimuons" in 
the 3//+track events (two entries per event) used as a control region for the analysis of events in 
the quadmuon category {R\). The distribution is compared with the prediction obtained from 
the full 2D background template, fitted for overall normalization only, (d): The invariant mass 
of the four "muons" in the 3f/+track control region for the quadmuon topology R\ compared 
with the prediction obtained from the data in the background-enriched region {2ji + 2tracks). 

p'^'' < 80 GeV/c (control region). The first sub-region is used to obtain a parametrization of 
the shape of the dimuon invariant mass distribution in background events. The data were pa- 
rameterized using a combination of Bernstein polynomials [16J and Crystall Ball functions flTl 
describing resonances, and the results are shown in Fig. [ija). To validate the template with 
data, we fit its shape to the observed data in the region 60 < pj^^ < 80 GeV/ c, allowing only 
the overall normalization to float in the fit. To account for the evolution of the resolution of 
mass measurement with pj^, an additional term was added to the uncertainties in the widths 
of resonances, and we used simulation to verify that the shape of the bulk of the background 
distribution is only a weak function of Pj^. The comparison shows good agreement as illus- 
trated in Fig.[2ta). The same template (with the uncertainty of mass resolution evolved to cover 
even higher p^^') is used to predict the shape of background events of topology R\ in the signal 
region pj'* > 80 GeV/c. 

The SM backgrounds in the two-dimuon topology R22 dominated by bb events with each b 
quark yielding a pair of muons. Because each b jet fragments independently, the background 
distribution in the {nii,m2) space of the two dimuon masses is a Cartesian product of the ID 
dimuon mass distribution with itself. However, because one of the dimuons contains the pj > 
15 GeV/ c muon that triggered the event, its dimuon mass spectrum is different from that of the 
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other dimuon. To account for this effect, we separately measure the shapes for the "trigger " and 
"non-trigger" dimuons. To model the trigger dimuon shape, we use single-dimuon events with 
further selections suppressing the non-bb backgrounds, fit to a parameterized functional form 
and subtract the residual contamination from Drell-Yan background, which does not contribute 
to the topology R22 (both subtracted and non-subtracted curves are shown in Fig. [T|b)). To 
match the kinematics of the two-dimuon events being modeled, the "other dimuon" shape is 
obtained using three-muon events with a dimuon recoiling off a trigger-quality muon and is 
shown in Fig. [l|c). To account properly for a contribution with both dimuons containing a 
trigger-quality muon in the final 2D template, an additional reweighing is applied in taking 
the Cartesian product of the two distributions. The template is validated using finally selected 
two-dimuon events in the off-diagonal part of the {mi,m2) distribution. Figure |2|b) shows a 
comparison of the invariant mass distribution of the dimuons in these events (note that there 
are two histogram entries per event), compared to the prediction derived from the template 
and fit to the data for overall normalization only. 

The quadmuon topology R\ has a small background contamination in which a b quark pro- 
duces two genuine muons and additional two muons are produced from non-muon tracks 
incorrectly matched to some of the genuine muon stubs. When identifying the two funda- 
mental dimuons within the group of four muons, both {ji, }i) + {trk, trk) and trk) + {ji, trk) 
pairings can occur, each having its own distinct 2D shape in the [m-[,m2) space. To model 
these events, we use single dimuon events and construct "pseudo muon jets" using two recon- 
structed muons and two non-muon tracks playing the role of misidentified muons. Selected 
events are separated into two subsets according to the type of pairing, each producing a 2D 
distribution for the invariant masses of the two pairs in the event. Figures [ifd), (e) and (f) 
show ID invariant mass distributions for (ji, {trk, trk), and trk)-ty-pe dimuons obtained 
from projections of the 2D distributions for the two types of events. In the R\ signal events, 
the identities of }i and trk are unknown, so the mass templates and the signal events are both 
symmetrized by randomly assigning dimuon masses to the horizontal and vertical axes of the 
2D distribution. The template is validated using a control region with three nearby muon 
candidates (K3), one of which is likely a misidentified hadron, and adding a non-muon track 
to play the role of a second misidentified muon. Figure |2]^c) compares the distribution of all 
"dimuons" in the 3/^(+track control sample (note two entries per event) compared to the pre- 
diction based on the full 2D template fitted to data for overall normalization only. Figure |2jd) 
makes a similar comparison but for the quadmuon invariant mass. Templates for higher order 
topologies are derived as combinations of the above methods. In all cases, the full posterior 
density functions for fit parameters including correlations were used in the final fit to account 
for the uncertainties in the background templates. 

3.4 Efficiencies and Systematic Uncertainties 

The shape of the invariant mass distribution for possible signal events was studied by com- 
paring the properties of events with dimuons from co, (p, ]/ip, and ^' resonances in data with 
the simulation predictions and extrapolating between the resonance masses. Because of the 
excellent resolution of the CMS tracker, signal shapes have narrow widths scaling with the 
mass of the resonances, with a slight dependence on the dimuon transverse momentum. For 
final fits, signal shape is parameterized using a Crystal Ball function with core resolution of 
a{m^^^) = 0.026 GeV/c^ + 0.0065 m,,,, for barrel (|?/W'| < 0.9) and 0.026 GeV/c^ + 0m3m^,fi for 
endcap (0.9 < < 2.4). To account for the unknown momentum range of the hypotheti- 

cal signal, the width of the Gaussian core was allowed to vary within 30%. Multidimensional 
distributions were obtained by taking appropriate Cartesian products. The uncertainties on 
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Figure 3: The ID and 2D invariant mass distributions of muon pairs for events in each signal 
region, compared with the expected background, (a): Events in the single dimuon topology 
(b): 10 events in the two-dimuon topology R22- ('^)- ^^'^ single "quadmuon" event in topology 
R4. (d): The invariant mass of all four muons for the same event. None of the events in the 
multi-dimensional topologies fall into the corridor along the diagonal (shown as dashed lines), 
which would indicate the presence of signal. The last plot is relevant for the special scenario 
with a cascade decay fl2 — > flifli with fn(a2) < 1m{a\), leading to the off-shell production of a\. 



the parameters of the function are obtained by quantifying the level of agreement between 
data and simulation, which is dominated by the statistical uncertainties. For dimuons with 
p'^'' < 150 GeV/ c, the reconstruction efficiency in the barrel region is nearly flat as a function 
of //, has an average value of 95 ± 1%, and is driven by the efficiency in finding and matching 
stubs in the muon system. The efficiency decreases to about 92% for m^^ close to 2m^ because of 
muon trajectories becoming nearly collinear. In the endcap region, there is a slight lowering of 
efficiency towards high | r] \ because of muon trajectories overlapping in the muon system. The 
systematic uncertainty on the efficiency in the endcap region is 3% as obtained from compar- 
isons of reconstruction efficiencies for simulated events with overlapping and non-overlapping 
muons, both with each other and with the data. For higher momentum dimuons, tracking ef- 
fects become important as the overlaps of the trajectories lead to high sharing of hits and de- 
creased efficiency. This effect is especially pronounced for m^,^, of 0.4—0.6 GeV/ c^, where the 
efficiency decreases to 85 ± 5% at ^ 250 GeV/c and to 75 ± 10% at p'^'' ~ 350 GeV/c be- 
cause of events in which muon trajectories bend towards each other in the magnetic field and 
remain close to each other for a substantial part of their paths through the tracker. For muon 
pairs with mj,^, outside of this range, the efficiency is nearly flat until p^'' « 250 GeV/ c, where it 
starts a slow descent as p^^ increases; it is about 85 ± 5% at ~ 350 GeV/ c. The uncertainties 
quoted account for possible inaccuracies in modelling the size of the tracker hit clusters in the 
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silicon tracker potentially affecting the rate of merging of hits from nearby tracks, which is the 
dominant cause of the inefficiency of track reconstruction. 

Reconstruction of quadmuons suffers more significantly from reconstruction failures in the 
muon system. In addition to the 3% inefficiency per muon, the probability to reconstruct all 
four muons in a quadmuon has a significant additional term related to small uninstrumented 
gaps between the wheels in the central part of the muon system. With an inefficiency of 8- 
10% per muon crossing the gap region and a significant probability for one or more muons to 
cross it, the average reconstruction efficiency for a quadmuon is 83 ± 3% {p^^^^ ^ 150 GeV/ c). 
For high-momentum quadmuons, tracking effects become more significant. In topologies with 
on-shell dimuons relevant to this analysis, the quadmuon efficiency is determined by the prob- 
ability to reconstruct separately each of the two low-mass dimuons comprising the quadmuon, 
as significant overlaps of trajectories for muons from different dimuons are rare. The momen- 
tum dependence of the quadmuon efficiency closely follows the reconstruction efficiency of a 
dimuon with momentum equal to the higher-momentum dimuon within the quadmuon. The 
reconstruction efficiency also has the same dependence on the invariant mass of the dimuons 
within the quadmuon. For a quadmuon with constituent dimuon mass of about 0.5 GeV/ 
(the worst-case scenario), the average efficiency at pj^^^^ ~ 250 GeV/ c is about 74% and has 
a systematic uncertainty of 2-3% due to the tracking efficiency. Higher-multiplicity muon jets 
have larger inefficiencies because of muon system reconstruction failures, but are reconstructed 
as lower-multiplicity jets. In the context of a specific model, the migration of events between 
the high-multiplicity topologies does not reduce the overall acceptance. Higher multiplicity 
muon jet reconstruction is less affected by tracking inefficiencies as the average momentum of 
dimuons is moderate, even for very high momentum (~ 400 GeV/ c) muon jets. 



4 Results 

The data in the regions used to search for signal (the "diagonal" regions of multi-dimensional 
distributions and the single-dimuon events with pj^^ > 80 GeV/ c in topology R\) were not 
looked at until all analysis selections and signal extraction techniques were finalized. When the 
signal regions were examined, no evidence for new resonance production was found within the 
sensitivity of this analysis. Figure |3] shows the observed data for select topologies and the ex- 
pected SM background contributions, which were obtained using the templates for each topol- 
ogy. The templates were normalized to the data in the off-diagonal regions in high-multiplicity 
topologies (all but R\) and directly fitted in a combined signal-plus-background fit for the 
single-dimuon case {R\). For topology J^^^/ 10 events were observed in the non-diagonal re- 
gion of the distribution with one event falling outside the m^/^, > 5 GeV/ c range and therefore 
not seen in Figure |3|b). The single off-diagonal event in the quadmuon topology R\ can be 
interpreted in the context of models with new light bosons produced off -shell via a process of 
the type fl2 ^i^i — ^ 4^. No events were found in any topologies with more than four muons. 

To interpret the results in a model-independent fashion, we set the 95% confidence level (CL) 
upper limits on the allowed production rate of the single-dimuon+X, quadmuon+X and two- 
dimuons+X topologies. To simplify further phenomenological interpretations of these results, 
the rate is defined as the cross section times appropriate branching fractions to produce a par- 
ticular signature times kinematic and geometrical acceptance, assuming an ideal detector. The 
"ideal detector acceptance" ttideai is the probability for an event to satisfy analysis selections 
for a specific topology taking into account momentum thresholds and f] ranges of muons, but 
ignoring all detector and reconstruction algorithm inefficiencies. For a given model, ^ideal can 
be evaluated using one of the standard event generators and requires no knowledge of the in- 
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strumental inefficiencies, which are absorbed into the experimental limit calculation.The upper 
limits are shown in Fig. Qa) and the main acceptance cuts, defining the ideal detector accep- 
tance for the three topologies, are given in the legend of the same figure. 

Limits are set using a Bayesian technique including integration over the systematic uncertain- 
ties in the signal and background shape parameterizations and the background normalization, 
which are treated as nuisance parameters using the true posterior density for the background 
normalization and log-normal priors in other cases. The limits shown in Fig. |4ja) account for 
systematic uncertainties in the knowledge of instrumental inefficiencies as well as the varia- 
tions of the reconstruction efficiencies for dimuon masses ranging from 0.25 to 5 GeV/c^ and 
the means of the muon jet (dimuon or quadmuon) momentum distributions up to 250 GeV/ c. 
These variations were treated as an additional systematic uncertainty and were 7% for the 
quadmuon topology, 20% for dimuon, and 35% for the two-dimuon topology. For higher mo- 
mentum muon jets, the limits become weaker because of diminishing track reconstruction effi- 
ciency. Other systematic uncertainties account for the precision in the luminosity measurement 
(4%), uncertainties in the efficiency of reconstructing and matching muon stubs to the tracks 
and triggering (1-4% depending on the topology), and track reconstruction efficiency (5-10%). 

The limits presented can be used as conservative bounds on higher multiplicity topologies, e.g. 
the limit for the quadmuon+X topology can be used to constrain the production rate of signals 
with muon jets containing four or more muons in it. These limits are conservative in assuming 
that reconstruction failures always remove events from consideration, whereas in reality events 
in which there are reconstruction failures usually enter another, lower multiplicity, topological 
category. The limit for the quadmuon+X topology can be reinterpreted for models predicting 
fl2 ^1^1 — ^ 4^ cascades with one or both fli bosons off-shell (if m(fl2) < 2m(fli)). In this case, 
the 95% CL upper limit on the production rate is the same as the quadmuon+X limit in Fig. Qa) 
except near m (02) ~ 3.5 GeV/ c^, the invariant mass of all four muons in the single off-diagonal 
event observed in topology R|. In that region, the limit can be conservatively estimated by 
multiplying the quadmuon+X limit by factor of 4.74/3.00 (the ratio of Bayesian 95% CL upper 
limits for the mean of a Poisson variable given one and zero observed events in the absence of 
backgrounds). 

To set limits for representative benchmark scenarios, we use two models of SUSY with the 
dark sector. The first SUSY model [7J assumes standard MSSM squark/ gluino production 
and cascades followed by the decay of the MSSM LSP into the dark sector Xi ^ ^dark/l^dark 
or Xi ~^ 7dark/(^dark/ where ;\;dark is the new cold dark matter candidate. The lightest MSSM 
neutralino mass is set to 400 GeV/c^ and the dark sector masses used are m(?Zdark) = 3, 
m(7dark) = 0-5, and m(;tdark) = 300 GeV/c^. The resulting limits are shown as a function of 
gluino mass m(g) {m{q) = m{g) /1.2 for the first two generations) in Fig.|4|b) for three different 
choices of branching fractions B(7daj.k ~^ VV)- The systematic uncertainties used are the same 
as for model-independent limits with two exceptions. First, we exclude uncertainties used to 
compensate for variations in efficiency with p^j because in this case the momentum spectra 
are known. Second, we include an uncertainty of 3% in the acceptance to account for uncer- 
tainties in proton parton distribution functions (PDF) by varying parameterizations within the 
CTEQ6.6 [18] family, and comparing the central values of CTEQ6.6L with NNPDF2.0 [19], and 
MSTW2008 [20J sets. As a reference, typical signal acceptance (the probability for an event to 
be reconstructed in one of the signal topologies) is of the order of 70-75% for m{q) =500-1000 
GeV/c-^ assuming B(7dark — > V-V) = 100%. For B(7dark — > MM) = ^0°/° ^^id 33%, the corre- 
sponding ranges change to 30-40% and 20-25%, respectively. 

The second model HI assumes squarks to be the MSSM LSP {m{q) = m{g)/1.2). Following 
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Figure 4: (a): 95% CL upper limits on the rate of the signals of new physics with leptonic 
jets for three topologies expressed as cross section times appropriate branching fractions to 
produce a particular signature (including branching fraction for 7dark ~^ V-j^) times kinematic 
and geometrical acceptance (muon momentum thresholds and rj ranges), assuming an ideal 
detector, (b): Limits for the Dark SUSY model [7J with the MSSM LSP decaying via x\ -> 
;^dark7dark + ;^dark^dark(^ TdarkTdark)/ with the ;^dark being the new LSP. (c) and (d): Limits on 
the model [8J where squark is the MSSM LSP decaying into a quark and a light hidden sector 
fermion decaying to a lighter hidden sector fermion with emission of either a dark photon (c) 
or a light dark-Higgs (d) decaying to two dark photons. 



production, squarks decay via q —> qn2, where n2 is a heavier dark sector fermion with the de- 
cay modes dominated by either n2 — > "iTdark i^i is a lighter dark fermion) or n2 — > Mi/2dark(~^ 
TdarkTdark) • For each of the two sub-models the limits on the production cross section are shown 
in Figs.|4](c) and (d) for three different choices of branching fractions B(7dark — ^ HJ^)- The dark 
sector masses are set to m(/2dark) = 1-2, m{hdark) = 0.5, m{n2) = 2, m(ni) = 0.5 GeV/c^. The 
cross section curves shown in Figs. |4]^c-d) assume universality of squark masses across three 
squark generations. The cross section is reduced if squark masses are not universal. The lim- 
its presented are the most stringent to date for models with dark SUSY sector from collider 
experiments. 



5 Summary 

A topology-based search for groups of collimated muons (muon jets) using a data sample cor- 
responding to an integrated luminosity of 35 pb^^ revealed no signal of new physics within 
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the sensitivity of the measurement. No events consistent with two or more decays of the 
new light boson species to a pair of muons were found in the data, and no excess over the 
SM backgrounds was observed for production of single high-p^ bosons decaying to pair of 
muons. Relaxing the assumption that the light bosons are on-shell and searching for possible 
fl2 — > flifli 4f/ cascades revealed one event, consistent with the background expectation. 

With these observations, we set limits on the production of new low-mass states decaying to 
pairs of muons that can be applied to a broad class of models predicting the leptonic jet signa- 
tures. These 95% CL limits exclude production of new physics in several event topologies with 
cr X B X ttideai in the range of 0.1-0.5 pb, where ftideal refers to the kinematic and geometrical 
acceptance of the analysis for an ideal detector. We also set model-dependent limits on several 
benchmark models predicting production of the new light states in the context of dark SUSY 
significantly extending the sensitivity of previous searches performed at the Tevatron. 
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